Conserved mass aggregation model with mass-dependent fragmentation 
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We study a conserved mass aggregation model with mass-dependent fragmentation in one dimen- 
sion. In the model, the whole mass m of a site isotropically diffuse with unit rate. With rate uj, a 
mass m"^ is fragmented from the site and moves to a randomly selected nearest neighbor site. Since 
the fragmented mass is smaller than the whole mass m of a site for A < 1, the on-site attractive 
interaction exists for the case. For A = 0, the model is known to undergo the condensation phase 
transitions from a fluid phase into a condensed phase as the density of total masses (p) increases 
beyond a critical density pc. For < A < 1, we numerically confirm for several values of lu that pc 
diverges with the system size L. Hence in thermodynamic limit, the condensed phase disappears 
and no transitions take place in one dimension. We also explain that there are no transitions in any 
dimensions. 
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Nonequilibrium phase transitions have been studied 
extensively during last decades and observed in various 
systems [l], 0, [3, As in equilibrium, nonequilibrium 
transitions can be classified into universality classes ac- 
cording to conservation laws and symmetries, which are 
characterized by several critical exponents. 

In this paper, we study condensation phase transi- 
tions. Nonequilibrium condensation phase transitions 
from fluid phase into condensed phase have been observed 
in a variety of phenomena ranging from traffic flow to 
polymer gels 1,1, 0, i, i, 0, [lE M, [H, [H, IH . In the 
steady state, a finite fraction of total particles condenses 
on a single site in condensed phase when the total parti- 
cle density p is increased beyond a certain critical value 
Pc- In fluid phase below pc, the particle number of each 
site fluctuates around p without the condensation ITU . 

The simplest model exhibiting the condensation is 
zero-range process (ZRP) [5]. In ZRP, many identical 
particles occupy sites on a lattice. Each site may contain 
an integer number of particles and one of these particles 
can hop to one of the nearest neighboring sites with a 
rate that depends on the number of particles on the site 
of departure. The chipping (single-particle dissociation) 
and aggregation processes of ZRP describe various con- 
densations such as traffic jam 0, bunching of buses 0, 
coalescence of shaken steel balls Q and the condensation 
on complex networks 0]. 

Another important class of condensation transitions 
emerges when the diffusion of the whole particles of a 
site is allowed in addition to the chipping and aggrega- 
tion. These processes arise in a variety of phenomena 
such as polymer gels [l^ , the formation of colloidal sus- 
pensions [ll| , river networks [H, [l^ and cloud formation 
[14| . Recently studied symmetric conserved-mass aggre- 
gation(SCA) model is the simplest one incorporating dif- 
fusion, chipping and aggregation upon contact [l^, [lB| • 



In one dimensional SCA model, the mass of site i 
moves either to site i — 1 or to site i + I with the unit 
rate, and then — > and mi±i — > mi±i + rrii. With 
the rate w, unit mass chips off from site i and moves to 
one of the nearest neighboring sites; mi ^ rm — 1 and 
mi±i — *■ TOiii -|- 1. As total masses are conserved, the 
conserved density p and to determine the phase of SCA 
model. The condensation transition arises from the com- 
petition between diffusion and chipping processes. The 
diffusion of masses tends to produce massive aggregates 
and consequently creates more vacant sites. The chip- 
ping of unit mass tends to prevent the formation of ag- 
gregates, so that it leads to a replenishment of the lower 
end of mass distribution. 

The single site mass distribution P(to), i.e., the prob- 
ability that a site has mass ni in the steady state, was 
shown to undergo phase transitions on regular lattices 
For a fixed lo, as p is varied across the critical value 
pMv), the behavior of P{m) for large m was found to be 



P{m) 



e -'"/'"* /5<Pc(w) 
m-^ p = pc{uj) (1) 

TO^"^ + infinite aggregate p > Pc{i^), 
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where pc is given as Pc(w) = \/uj + 1 — 1. pc and r 
are^shown to be independent of the spatial dimension 
d [16]. The tail of the mass distribution changes from 
exponential to an algebraic decay as p approaches pc 
from below. As one further increases p beyond pc, this 
asymptotic algebraic part of the critical distribution re- 
mains unchanged, but in addition an infinite aggregate 
forms. This means that all the additional mass {p— pc)L'^ 
condenses onto a single site and does not disturb the 
background critical distribution. The a; = oo case corre- 
sponds to ZRP with a constant chipping rate, and then 
there is no condensation transitions on regular lattices 
0. The critical exponent r is same everywhere on the 
critical line Pc(^)- Recent studies showed that the ex- 
istence of the condensation transitions depends on the 
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spatial disorder , the symmetry of moving directions 
[l^ , the constraint of diffusion rate [T^ , and the under- 
lying network structure (20l |. 

In more general situations, as pointed out in Ref. [TH|, 
the diffusion rate of mass and the chipped mass will 
depend on the mass of a departure site. When the 
diffusion rate depends on mass as D{m) ~ m~" with 
a > [l^, the condensed phase disappeared in the 
thermodynamic limit. However finite systems undergo 
the condensation transitions at a certain critical density 
Pc which algebraically diverges with the system size L. 
On the other hand, the mass-dependent fragmentation 
or chipping is also important in general physical situa- 
tion. For instance, in gelation phenomena, there is no 
reason that only one monomer at the end of a poly- 
mer is chipped regardless of the mass of polymers. It 
is very likely that more monomers can be chipped for 
more massive polymers. Hence it is natural to study the 
effect of mass-dependent fragmentation on the conden- 
sation phase transitions. In this paper, we investigate 
the condensation transitions of SCA models in which the 
fragmented mass in a chipping process depends on the 
whole mass of the departure site. 

We consider SCA model with mass-dependent frag- 
mentation in one dimension. In the model, the frag- 
mented mass in the chipping process is given as Sm — 
Am^ . By chipping process with the rate oj, the mass 
Sm randomly moves to the one of the nearest neighbors, 
rrii rrii — Srrii and mi±i mi±i + Srrii. With the unit 
rate, mass rrii randomly diffuses to the one of the nearest 
neighbors, rrii and rni±i — > mi±i + rrii. 

Since the chipped mass is smaller than the whole mass 
of a site for A < 1, the attractive interaction exists be- 
tween masses on the same site. The A = case is just 
SCA model of Ref. [l^. The A = 1 is a trivial point, but 
two extreme situations emerge according to the values of 
A. For yl = or 1, masses only diffuse without chipping 
so the complete condensation takes place. On the other 
hand, for < A < 1, the fragmented mass is proportional 
to the whole mass, which means that any aggregates can- 
not exit. Hence, the system is always in the fluid phase. 
The A = 1 and < A < 1 case corresponds to the ran- 
dom fragmentation model in which P{m) always expo- 
nentially decays [2l[. On the other hand, for A < 0, the 
complete condensation should always occur because the 
time needed to dissipate an aggregate, if exists, exponen- 
tially increases with the mass of an aggregate. Therefore, 
we focus on the range of < A < 1. 

We perform A^Ionte Carlo simulations on a ring of size 
L. With random initial distributions of total mass pL, 
the system is allowed to reach the steady states. We run 
simulations typically up to t = 10® time steps for size L 
up to 2048. In the steady states, we measure the single- 
site mass distribution P{m) for various values of A and to. 
Fig.l shows Pirn) of A = 0.1 and 0.25 with p = 20 and 
uj = 0.1 for L up to 2048. For both A values, an aggregate 
with mass (p — Pc)L forms and its probability decreases 
as 1/L. Pc is a certain critical density above which an 
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FIG. 1: The double logarithmic plots of P(m) for A = 0.1 
(a) and 0.25 (b) with u — 0.1. Insets shows the scaling plot 
m-^P{m) with r = 1.84 (a) and 1.62 (b) for the data of 
L — 2048. In each panel, the solid line corresponds to p = 0.01 
with L = 2048. Symbols are the data of p = 20 with L up to 
2048. 



aggregate forms. P{m) shows that in finite systems, the 
condensation phase transitions occurs at finite pc for A > 
0. For other values of uj, the condensation transitions 
also occur in finite systems, pc increases with lu while r 
is the same for all w. 

If there exists the condensation transitions in the ther- 
modynamic limit, r should satisfy r > 2. We measure 
the exponent r from the scaling plot of m'^P(m). For 
L = 2048 with w = 0.1, we estimate r as 1.84(10) for 
A = 0.1 and 1.62(10) for A 0.25. For other values of cj, 
we obtain the similar r within errors. However, the finite 
size effect appears at very small mass even for large L, for 
example, at about m = 100 for L = 2048. Even though 
the precise measurement of r is difficult due to the strong 
finite size effect, our estimates of t are smaller than 2. It 
means that the condensation transitions take place only 
in the finite systems and disappear in the thermodynamic 
limit. In what follows, we discuss that no condensation 
transitions occurs in the thermodynamic limit. 

In infinite systems with finite density ( or L — > cxd ), 
T must be greater than 2 since the first moment < m > 
is the density. Otherwise, the first moment diverges as 
InL for T = 2 and for r < 2 respectively. Hence 
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FIG, 2: The double logarithmic plots of pc{L) for A — 0,1 
(a) and 0,25 (b) with oj = 0,1 and p — 20, The slope of line 
in each panel is /3 = 0,1 (a) and 0,25 (b). 



T < 2 for A > contradict with the conservation of total 
density. Similar finite size effect with t < 2 has been 
observed in asymmetric CA [3] and SCA with mass- 
dependent diffusion rate [l^ , pc diverges with L as pc ^ 
InL for asymmetric CA and for the mass-dependent 
diffusion rate respectively. In the thermodynamic limit 
or L — !■ cx), no transitions occurs due to pc = oo and a 
system is always in the fluid phase. As a result, r < 2 
for A > implies that pc = oo and the condensation 
transitions disappear in the thermodynamic limit. Only 
finite systems undergo the condensation transitions at a 
finite Pc- 

To find the size dependence of Pc, we measure pc for 
each L as follows. In the steady states, P{m) scales as 

M 



P{m) = m-^J{m/L'^) + -5{m - {p - p,)L) , (2) 

where is the crossover exponent [iGj . From the conser- 
vation of total masses, pc is given as pc = P — Poo, where 
Poo is the density of an aggregate. From the fact that 
the background distribution does not change for p> Pc, 
one can estimate Pc via " mP{m)dm in the condensed 
phase, where toq is the cut-off mass at which the back- 
ground distribution terminates. Using this method, we 



confirmed that our estimate for pc of the normal SCA 
model 15] converges to the the exact pc as L increases 
(not shown). In Fig. 2, we plot pc as a function of L 
for A = 0.1 and 0.25 with lu — 0.1, where we measure pc 
using P{m) for p ~ 20. Insets show the local slope Peff 
defined as Peff = \n[p{L)/ p{L/b)]/ In b with 6 = 2. As 
one can see, pc diverges with L as with (3 = 0.11(3) 
for A = 0.1 and /3 = 0.25(5) for A = 0.25. Since the frag- 
mented masses increases for A — > 1, pc should diverge 
faster for larger A. Hence we conclude that there are no 
condensation transitions in the thermodynamic limit for 
< A < 1. 

The scaling relation between three exponents, r, (j) and 
/? is obtained by calculating pc = J™" mP{m)dm for 
P^ Pc where ttIc is the crossover mass over which P[m) 
decays exponentially due to the finite size effect. The 
lower limit is the mass from which P(m) begins to decay 
algebraically to the upper limit TOc. We set the lower 
limit of the integral to unit for convenience. Since rric 
scales with L as rric ^ L't', Pc scales as L'^(2-r)^ Hence 
f3 is given as P = (p{2 — r). From the normalization 
condition of _P(m), one obtains another scaling relation 
0(t— 1) = 1 Jl6f|. Hence we have the two scaling relations, 



P^H2- t) 



1) = 1. 



(3) 



From the two relations, r is given as r = (/3-|-2)/(/3+ 1). 
With the estimates of /3, we obtain r = 1.91 for A — 
0.1 and 1.8 for A = 0.25 respectively. Since the direct 
measurement of t crucially depends on the choice of the 
scaling region, the estimate of t is very subjective. On 
the other hand, pc exhibits the nice scaling with L as 
shown in Fig. 2, the prediction of r by the scaling relation 
should be more precise. Hence we predict r — 1.91(3) for 
A = 0.1 and 1.80(3) for A = 0.25. 

Recent studies showed that SCA models is well de- 
scribed by mean field theory [H [H, [13, [S El • 

Hence 

we expect that the present model is also well described 
by the mean field theory. To prove no transitions in the 
thermodynamic limit, we first assume the existence of the 
condensed phase, so that there is an aggregate. We also 
assume that all sites except the site occupied by an aggre- 
gate are occupied by mass Pc, which reflects the uniform 
background distribution with mass pc- Then, for the ex- 
istence of an aggregate in the steady states, the mass loss 
of an aggregate by chipping should be equal to the gain 
from the background. However the mass of an aggregate 
diverges with L, so the chipped mass of an aggregate also 
diverges as L'^. On the other hand, the background dis- 
tribution should not change in the condensed phase, the 
gain from the background is finite and independent of L. 
As a result, the loss diverges in the thermodynamic limit 
while the gain is finite. Therefore an aggregate cannot 
exist in the thermodynamic limit, which contradicts with 
the assumption of the existence of an aggregate. So the 
condensed phase cannot exist and the system is always in 
the fluid phase in the limit L ^ oo. Since the argument 
is based on the mean field, we expect no transitions in 



4 



any dimensions. 

In summary, we investigate the condensation phase 
transitions of SCA model with a mass-dependent frag- 
mentation in one dimension. We numericaUy find that 
condensation transitions take place only in finite sys- 
tems. The critical density diverges with system size L 
as because of r < 2 for < A < 1. Hence in the 
thermodynamic limit of L — s- cx), no condensation transi- 
tions occur, so that a system is always in the fluid phase 



for A > 0. Based on the mean field argument, we con- 
clude that there are no condensation transitions in any 
dimensions. 
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